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ABSTRACT 

We study the Local Group spiral galaxy M33 to investigate how the observed scaling between the 
(kpc-averaged) surface density of molecular gas (£h2) and recent star formation rate (Ssfr) relates 
to individual star-forming regions. To do this, we measure the ratio of CO emission to extinction- 
corrected Ha emission in apertures of varying sizes centered both on peaks of CO and Ha emission. 
We parameterize this ratio as a molecular gas (H 2 ) depletion time (rd cp )- On large (kpc) scales, 
our results are consistent with a molecular star formation law (Xsfr ~ ^112) with b ~ 1.1 — 1.5 
and a median 7"d ep ~ 1 Gyr, with no dependence on type of region targeted. Below these scales, 
Td C p is a strong function of adopted angular scale and the type of region that is targeted. Small 
(< 300 pc) apertures centered on CO peaks have very long Td op (i.e., high CO-to-Ha flux ratio) and 
small apertures targeted toward Ha peaks have very short Td ep - This implies that the star formation 
law observed on kpc scales breaks down once one reaches aperture sizes of < 300 pc. For our smallest 
apertures (75 pc), the difference in Td ep between the two types of regions is more than one order of 
magnitude. This scale behavior emerges from averaging over star-forming regions with a wide range 
of CO-to-Ha ratios with the natural consquence that the breakdown in the star formation law is a 
function of the surface density of the regions studied. We consider the evolution of individual regions 
the most likely driver for region-to-region differences in Td op (and thus the CO-to-Ha ratio). 

Subject headings: Galaxies: individual (M33) — Galaxies: ISM — H II regions — ISM: clouds - 
Stars: formation 



1. INTRODUCTION 

The observed correlation between gas and star forma- 
tion rate surface densities (the 'star formation law') is 
one of the most widely used scaling relations in extra - 
galactic astronomy (e.g., [Schmidt 1959i lR"ennicu 
However, its connection to the fundamental units of star 
formation, molecular clouds and young stellar clusters, 
remains poorly understood. On the one hand, averaged 
over substantial areas of a galaxy, the surface density of 
gas correla tes well with th e amount of recently formed 
stars (e.g.. !Kennicuttiri998| ). On the other hand, in the 
Milky Way giant molecular clouds (GMCs), the birth- 
place of most stars, and H II regions, the ionized ISM 
regions around (young) massive stars, are observed to 
be distinct objects. While they are often found near 
one another, the radiation fields, stellar winds, and ulti- 
mately supernovae make H II regions and young clus- 
ters hostile to their parent clouds on small (~10 pc) 
scales. Thus, while correlated on galactic scales, young 
stars and molecular gas are in fact anti-correlated on 
very small scales. The details of the transition between 
these two regimes remain largely unexplored (though see 
lEvans et al.ll2009t) . 

Recent observations of nearby galaxies have identi- 
fied a particularly tight correlation between the distri- 
butions of mole cular gas (H2) and recent star formation 
on -kpc scales (jMurgia et alj|2002t FWong fc Blitzll2002t 
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Kennicutt et al.l 120071: IBigiel et all 120081: ILerov et al.l 
2008!: IWilson" et al. l ^oMT While the exact details of the 
relation are still somewhat uncertain, in the disks of spi- 
ral galaxies the parametrization seems to be a power law, 
^sfr = a ^H2J w hh power law index b m 1 — 1.7 and co- 
efficient a corresponding to H2 depletion times of —2 Gyr 
in normal spirals. 

Both parts of this relation, the surface densities of H 2 
and recent star formation, resolve into discrete objects: 
GMCs and H II regions, young associations, and clusters. 
In this paper, we investigate whether the —kpc H2-SFR 
relation is a property of these individual regions or a con- 
sequence of averaging over large portions of a galactic 
disk (and the accompanying range of evolutionary states 
and physical properties). To do so, we compare CO and 
extinction-corrected Ha at high spatial resolution in the 
nearby spiral galaxy M33. We examine how the ratio of 
CO-to-Ha changes as a function of region targeted and 
spatial scale. M33 is a natural target for this experi- 
ment: it has favorable orientation and is close enough 
that peaks in the CO and Ha maps ap proximately cor- 
respo nd to individual massive GMCs (Rosolowsky et all 
2007]) and H II regions (jHodge et al.ll2002D . 

Perhaps not surprisingly, we find that the ratio of CO 
luminosity (a measure of the molecular gas mass) to 
extinction-corrected Ha flux (a measure of the star for- 
mation rate) depends on the choice of aperture and spa- 
tial scale of the observations. After describing how we 
estimate H2 masses and the recent star formation rate 
(Section [2]) and outlining our methodology (Section |3]), 
we show the dependence of the depletion times on spatial 
scale and region targeted (Section |4j . Then we explore 
physical explanations for these results (Section [5]). 
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2. DATA 

We require the distributions of H 2 and recently formed 
stars which we trace via CO emission and a combination 
of Ha and IR emission, respectively. 

2.1. Molecular Gas from CO Data 

Star- forming clouds consist mainly of H2, which can- 
not be directly observed under typical conditions. In- 
stead, H2 is usually traced via emission from the sec- 
ond most common molecule, CO. We follow this ap- 
proach, estimating H2 ma s ses fro m the CO J = 1 — 
data of iRosolowsky et al.l (12007|). which combines the 
BIMA (interferometric) data of fengargiola et al.l (120031) 
and th e FCRAO 14 m (single-dish) data of lHever et al.l 
(2004). The resolution of the merged data cube is 
13" x 2.03 km s _1 with a median la noise of 240 mK 
(~2.1 M fl pc~ 2 for our adopted Xco)- IRosolowsky et al.l 
(2007) showed that this c ombined cube recovers the flux 
of the lHever etall (12001 FCRAO data. 

We convert integrated CO intensities into molecu- 
lar gas surface densities assuming Aco = 2.0 x 10 20 



cm (K km s 



This is approximately the Milky 



W ay conversion facto r and a grees well with work on M33 
bv IRosolowsky et al.l ([2003D . For this X co : 



E H2 [M pc- 2 ] = 4.4 7 CO [K km s" 1 ] , 



(1) 



where Iqo is the integrated CO intensity over the line of 
sight and Sh2 is the mass surface density of molecular 
gas, including helium. 

The data cover a wide bandpass, only a small portion 
of which contains the CO line. As a result, direct inte- 
gration of the cube over all velocities produces an unnec- 
essarily noisy map. Therefore, we "mask" the data, iden- 
tifying the velocity range likely to contain the CO line 
along each line of sight. We integrate over all channels 
with ±25 km s -1 of the local mean H I ve locity (using 
the data from iDeul fc van der Hulsti[l~98l . To ensure 
that this does not miss any significant emission, we also 
convolve the original CO cube to 30" resolution and then 
identify all regions above 3cr in 2 consecutive channels. 
Any region within or near such a region is also included 
in the mask. We blank all parts of the data cube that 
do not meet either criteria and then integrate along the 
velocity axis to produce an integrated CO intensity map. 
Figure Q] shows this map at full resolution {middle left) 
and smoothed to ^45" resolution (top left) to increase the 
SNR and highlight extended emission. The noise in the 
integrated intensity map varies with position but typical 
la values are 8 — 10 M Q pc -2 ; the dynamic range (peak 
SNR) is ^20. The 3er mass sensitivity in an individual 
resolution element is ~10 5 M . 

2.2. Recent Star Formation from Ha and IR Data 

We trace the distribution of recent star formation us- 
ing Ha emission, which is driven by ionizing photons pro- 
duced almost exclusively in very young (massive) stars. 
We account for extinction by combining Ha and in- 
frared (24/i m) emission, a powe rful te chnique demon- 
strated by iCalzetti et al.l ([20071 ) and iKennicutt et al.l 
(|2007l I2009I ). Assuming continuous star formation over 
the past 100 Myr and studying a set of extragalactic star- 
forming regions, Calz etti et al.l ([20071 ) found the recent 
star formation rate (SFR) to be 



SFR [M yr" 1 ] = 5.3xl0" 42 [L (Ha) + 0.031 L (24/im)] , 

(2) 

where L(Ra) and L(24/im) = i/L„(24/im) are the lumi- 
nosities of a region in Ha emission and at 24/im, mea- 
sured in erg s . 

The assumption of continuous star formation is 
certainly inapplicable to individual regions, which 
are better described by instantaneous bursts (e.g., 
IRelano fc Kennicuttll2009l ). We only report averages of a 
large set (~150) of regions, which together constitute a 
large part of M33's total Ha, and argue that this justifies 
the application of Equation [5] (see Section |3~2"|) . In any 
case SFR units allow ready comparison to previous work. 

2.2.1. Ha Data 

We use t he n arrow-band Ha image obtained by 
IGreenawaltl (|1998l ) with the KPNO 0.6 m telescope. The 
reduction, continuum subt raction, and other details o f 
these data are described bv lHoopes fc Walterbosl {2000). 
Before combination with the IR map, we correct the 
Ha map for Ga lactic extinction usi ng a reddening of 
E(B-V) = 0.042 (ISchleeel et al.l [1991 and a ratio of Ha 
narrow band extinction to reddening of R(Ha) = 2.33. 

Studies of M33 and other nearby galaxies find typi- 
cally ~40% of t he Ha emission to come from "diffuse ion- 
ized gas" (DIG.IHoopes fc Walterbol[2000l: iThilker et al.l 
I2002T iHoopes fc Walterbosl 120031: iThilker et al.l 120051) . 
The origin of this emission is still debated; it may be 
powered by leaked photons from bright H II regions or 
it may arise "in situ" from isolated massive stars. We 
choose to remove this diffuse component from the Ha 
map and any discussion of Ha emission in the following 
analysis refers to this DIG-subtracted map (we assess 
the impact of this step in the Appendix). We do so us- 
ing the following method from IGreenawaltl ([1998). We 
begin by median filtering the Ha map with a 900 pc ker- 
nel. We then identify H II regions as areas in the original 
map that exceed the median-filtered map by an emission 
measure of 50 pc cm~ 6 (outlined by black contours in the 
top right panel of Figure [T]) . We blank these regions in 
the original map and smooth to get an estimate of DIG 
emission towards the H II regions. Our working map con- 
sists of only emission from the H II regions after the DIG 
foreground has been subtracted. The integrated Ha flux 
(inside R ga i = 4.5 kpc) allocated to the diffuse map is 
1.05 x 10 40 erg s -1 (44%) while the part allocated to the 
DIG-subtracted, H Il-region map is 1.35 x 10 40 erg s _1 
(56%), in good agreement with previous results on M33 
and other nearby galaxies. 

2.2.2. IR Data 

We measure IR intensities from 24 /xm maps obtained 
by the Spitzer Space Teles cope (PI: iGehrz et al.l 120051 
see also IVerlev et al.l |2007). The data were reduced 
by K. Gordon (200 9, private communication) following 
IGordon et al.l ([2005D . Spitzer' s point spread function at 
24 fim is ~6", well below our smallest aperture size 
(~18") and so is not a large concern. 

As with the Ha image, the 24 /j,m map includes a 
substantial fraction of diffuse emission — infrared cir- 
rus heated by an older population, emission from low- 
mass star-forming regions, and dust heated leakage from 
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Fig. 1. — CO and recent star formation in M33. Top left: Masked, integrated CO intensity smoothed to ~45" resolution to enhance the 
SNR. The black ellipse shows iJ ga i = 4.5 kpc (^0.6 r2s); we carry out our analysis inside this radius. Top right: Ho emission corrected for 
extinction via combination with mid-IR, emission. A black contour outlines regions that remain after the subtraction of the diffuse ionized 
gas (see text). Middle l eft: Positions of our CO peaks (blue diamonds) plotted on the full resolution CO map along with the cataloged 
positions of CMCs from Rosolowskv et al. (2007, green squares). Middle Right: Positions of our Ha peaks (blue diamonds) along with the 
150 most luminous H II regions cataloged by Hodge ct al. (2002, green squares). Bottom left: Relative distributions of bright CO (red) and 
Ha (blue) emission and the area of overlap (purple), inflated by 6" for display reasons. Bottom right: Histograms of distance from each 
peak to the nearest peak of the same type (red, i.e., Ha to Ha or CO to CO) and of the other type (blue, CO to Ha and vice versa). 
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nearby H II regions, with a mi nor contribution from pho- 
tospheric emission of old stars. iVerlev et all ([2007) argue 
that this diffuse emission accounts for ~2/3 of all 24 fim 
emission in M33. To isolate 24 fim. emission originating 
directly from H II regions, we follow a similar approach 
to that used to remove DIG from the Ha map. The key 
difference is that instead of trying to identify all 24 /im 
bright sources by filtering and applying a cut-off to the 
24 /im emission, we use the existing locations of H II re- 
gions to isolate any local 24 /im excess associated with 
H II regions. We extinction-correct the DIG-subtracted 
Ha emission using only this local excess in 24 /im emis- 
sion. The total integrated flux at 24 fim (R ga i < 4.5 kpc) 
is 3.92 x 10 41 erg s _1 , the fraction of DIG-subtraced 
24 /im inside the H II region mask is 1.63 x 10 41 erg s _1 
(42%). The the 24 /zm correction implies Ha extinctions 
of Ahq ~ 0.3 — 0.4 magnitudes. 

3. METHODOLOGY 

To quantify the scale-dependence of the molecular 
star formation law, we measure the H2 depletion timc0, 
Tdop = Eh2/Ssfr, for apertures centered on bright CO 
and Ha peaks. We treat the two types of peaks sepa- 
rately and vary the sizes of the apertures used. In this 
way we simulate a continuum of observations ranging 
from nearly an entire galaxy ( > 1 kpc apertures) to stud- 
ies of (almost) individual GMCs or H II regions (75 pc 
apertures). The CO data limit this analysis to galacto- 
centric radius < 4.5 kpc (~0.6 r 25 ). 

3.1. Identifying CO and Ha Peaks 

We employ a simple algorithm to identify bright re- 
gions in the DIG-subtracted, extinction-corrected Ha 
map and the integrated CO intensity map. This auto- 
mated approach allows us to use the same technique on 
both maps to find peaks matched in scale to our small- 
est aperture (75 pc). It is also easily reproducible and 
extensible to other galaxies. 

This algorithm operates as follows: We identify all con- 
tiguous regions above a certain intensity — the local 3cr 
in the CO map and ~1.9 x 10 40 erg s _1 kpc -2 in the cor- 
rected Ha map (~0.1 M© yr _1 kpc -2 following Equation 
[5]). We reject small regions (area less than ~110 arcsec 2 , 
which correspond to ^1800 pc 2 at the distance of M33) as 
potentially spurious; the remaining regions are expanded 
by 20" (^80 pc) in radius to include any low intensity 
envelopes. The positions on which we center our aper- 
tures are then the intensity- weighted average position of 
each distinct region. 

We find 172 CO regions and 154 Ha regions. Strictly 
speaking, these are discrete, significant emission features 
at ^50 pc resolution. At this resolution, there is a close 
but not perfect match between these peaks and the real 
physical structures in the two maps — GMCs and H II 
regions. Figure [1] (middle panels) sho ws our peaks along 
with the cataloged positi ons of GMCs (IR psolowsk vet al.l 
I2007D and H II regions (jHodge et al.l 12002ft . There is 
a good correspondence, with > 80% of the 149 known 
GMCs and the 150 brightest H II regions lying within 
~6" (3 pixels) of one of our regions. 

5 We emphasize that Tdep maps directly to observablcs. It is 
proportional to the ratio of CO to extinction-correct Ha emission. 



TABLE 1 

T dcp as Function of Peak and Scale 



Scale 

(PC) 


Depletion Time (Gyr) 
centered on CO centered on Ha 


<iV) a 


1200 


1.1 ± 0.1 


0.9 ±0.1 


16.2 


600 


1.2 ±0.3 


1.0 ±0.1 


5.2 


300 


2.5 ± 0.5 


0.64 ± 0.05 


2.1 


150 


4.9 ±0.9 


0.41 ±0.04 


1.4 


75 


8.6 ± 2.1 


0.25 ±0.02 


1.1 



a Typical number of individual CO or Ha peaks in- 
side an aperture. 



3.2. Measuring Depletion Times 

For a series of scales d, we center an aperture of diam- 
eter d on each CO and Ha peak and then measure fluxes 
within that aperture to obtain a mass of H 2 (M^) and a 
star formation rate (SFR). We then compute the median 
H 2 depletion time for the whole set of apertures. We do 
this for scales d = 1200, 600, 300, 150 and 75 pc and 
record results separately for apertures centered on CO 
and Ha peaks. 

At larger spatial scales, apertures centered on different 
peaks overlap (because the average spacing between CO 
and Ha peaks is less than the aperture size). To account 
for this, we measure only a subset of apertures chosen so 
that at least 80% of the selected area belongs only to one 
aperture targeting a given peak type (CO or Ha) at one 
time. 

While we center on particular peaks, we integrate over 
all emission in our maps within the aperture. At the 
smallest scales we probe (75 pc), this emission will arise 
mostly — but not exclusively — from the target region. 
At progressively larger scales, we will integrate over an 
increasing number of other regions. 

3.3. Uncertainties 

We estimate the uncertainty in our measurements us- 
ing a Monte-Carlo analysis. For the high-SNR Ha and 
24 /im maps, we add realistic noise maps to the observed 
"true" maps and repeat the identification and removal 
of DIG emission using smoothing kernels and emission 
measure cuts perturbed from the values in Section 12.2.11 
by ±25%. The low SNR of the CO data requires a more 
complex analysis. We assume that all regions with sur- 
face densities above 10 M Q pc -2 (~1.5 a) in the inte- 
grated CO map contain true signal. We generate a noise 
map correlated on the (13") spatial scale of our CO data 
and scale this noise map according to the square root of 
the number of channels along each line of sight in our 
masked CO cube (typically 5 — 7). Then we add all 
emission from the pixels above 10 M Q pc -2 . Finally, we 
re-identify peaks in the new maps and re-measure Mh 2 
and SFR in each region. We repeat this process 100 
times; the scatter in Td cp across these repetitions is our 
uncertainty estimate. 

4. RESULTS 

Figure [5] shows a well-known result for M33. There is a 
strong correlation (rank correlation coefficient of r « 0.8) 
between the surface densities of SFR and H 2 at 1200 pc 
scales. Power law fits to the different samples (types 
of peaks) and Monte Carlo iterations yield H 2 depletion 



Scale Dependence of Td ep in M33 



5 



o.i 



0.01 



0.0001 



o.i 



Apertures Centered on CO Peaks 
Apertures Centered on Ha Peaks 



. 4* . 



• # 

*■ • 
★* 



0.1 



1 10 
H 2 Surface Density [M pc~ 2 ] 



0.01 



0.001 



Apertures Centered on CO 
Apertures Centered on Ha 


Peaks •/' 
Peaks.* 




75^ 


O J" 




': i5o 


* 






300 + 








600 * ' 
/ 1200| 


t 

300 


f f 

150 / 

/ 75 




600 









i io ioo 

Median H 2 Surface Density [M pe~ 2 ] 



Fig. 2. — The relationship between SFR surface density (j/-axis) 
and H2 surface density (rr-axis) at large spatial scales. Red points 
show measurements integrated over (one independent subset of) 
1200 pc apertures centered on CO peaks. Blue points show similar 
measurements centered on Ha peaks. Dashed lines indicate fixed 
H2 depletion times ranging from 0.1 Gyr in the upper left to 10 Gyr 
in the lower right. 
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Fig. 3. — Scale dependence of the H2 depletion time, r^ cp , in 
M33. The y-axis shows the logarithm of the median H2 deple- 
tion time for apertures of different diameters (rr-axis) centered on 
CO peaks (red) and He* peaks (blue). Error bars correspond to 
uncertainty in the median estimated via a Monte-Carlo analysis. 
Dashed lines show expectation for simply averaging together two 
populations of regions in different states (Section [5j. 

times, Td C p = M112/SFR, of ~1 Gyr (at H 2 surface den- 
sities of 3 Mq pc -2 ) and power law indices of ~1.1 — 1.5. 
These results (modulo some renorm alization due to dif - 
fere nt assumptio n s) ma tch those of iHeyer et al.l (|2004l ) 
and iVerlev et al.l ((2010) in their studies of the star for- 
mation law in M33. The important point here is that 
there is good evidence for an internal H2-SFR surface 
density relation in M33. 

We plot the median Td ep as a function of scale (aperture 
size) in Figure |3l giving results for apertures centered on 



Fig. 4. — Scale dependence of the location of data in the star 
formation law parameter space. Red points show the median Sspr 
(y-axis) and £h2 (x-axis) for apertures centered on CO peaks. Blue 
stars show the same for apertures centered on Ha peaks. Dashed 
lines as in Figure [2] 



CO (red circles) and Ha (blue stars) peaks seperately. 
For the largest scales, we find a similar Td op for both sets 
of apertures (as was evident from Figure [5]). Going to 
smaller aperture sizes, Td ep becomes a strong function of 
scale and type of peak targeted. Small apertures centered 
on CO peaks have very long Td ep (up to 10 Gyr). Small 
apertures targeted toward Ha peaks have very short Td ep 
(0.3 Gyr). This may not be surprising, given the expec- 
tations that we outlined in Section [T] and the distinctness 
of the bright Ha and CO distributions seen in the lower 
left panel of Figure [I] but the dramatic difference as one 
goes from ~kpc to ~100 pc scales is nonetheless striking. 

A few caveats apply to Figure [3] First, in subtracting 
the diffuse emission (DIG) from the Ha map, we removed 
^40% of the flux. This could easily include faint regions 
associated with CO peaks, which instead show up as ze- 
ros in our map. Perhaps more importantly, we use the 
24 /mi map only to correct the DIG-subtraced Ha map 
for extinction. Any completely embedded star formation 
will therefore be missed. For both of these reasons, the 
SFR associated with the red points, while it represents 
our best guess, may be biased somewhat low and cer- 
tainly reflects emission from relatively evolved regions — 
those regions that have Ha fluxes above our DIG-cutoff 
value. There is no similar effect for the CO map. 

Figure [3] implies that there is substantial movement of 
points in the star formation law parameter space as we 
zoom in to higher resolution on one set of peaks or an- 
other. Figure 0] shows this behavior, plotting the median 
Ssfr and median Sh2 for each set of apertures (N.B., 
the ratio of median Sh2 to median Ssfr does not have 
to be identical to the median Td op ; the difference is usu- 
ally ^ 30%). We plot only medians because individual 
data are extremely uncertain, include many upper limits, 
and because we are primarily interested in the systematic 
effects of resolution on data in this parameter space. 

Apertures centered on CO peaks (red points) have 
approximately constant Esfr, regardless of resolution. 
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This can be explained if emission in the Ha map is ho- 
mogeneously distributed as compared to the position of 
CO peaks. Meanwhile there is a strong change in Eh2 
for decreasing aperture sizes on the same peaks; Eh2 
goes up as the bright region fills more and more of the 
aperture. A similar effect can be seen for the Ha (blue 
stars), though there is more evolution in Ejj2 with in- 
creasing resolution because most bright Ha regions also 
show some excess in CO emission. 

5. DISCUSSION 

Figure [3] shows that by zooming in on an individual 
star-forming region, one loses the ability to recover the 
star formation law observed on large scales. For aper- 
tures < 300 pc in size, the relative amounts of CO emis- 
sion and Ha intensity vary systematically as a function 
of scale and what type of region one focuses on. An- 
other simple way to put this, demonstrated in Figure^ is 
that scatter orthogonal to the SFR-H2 relation increases 
with increasing resolution. Eventually this washes out 
the scaling seen on large scales and the star formation 
law may be said to "break down" . 

What is the origin of this scale dependence? In prin- 
ciple one can imagine at least six sources of scale depen- 
dence in the star formation law: 

1. Statistical fluctuations due to noise in the maps. 

2. Feedback effects of stars on their parent clouds. 

3. Drift of young stars from their parent clouds. 

4. Region-to-region variations in the efficiency of star 
formation. 

5. Time-evolution of individual regions. 

6. Region-to-region variations in how observables map 
to physical quantities. 

Our observations are unlikely to be driven by any of 
the first three effects. In principle, statistical fluctua- 
tions could drive the identification of Ha and CO peaks 
leading to a signal similar to Figure [3] purely from noise. 
However, our Monte Carlo calculations, the overall SNR 
in the maps, and the match to previous region identifi- 
cations make it clear that this is not the case. 

Photoionization by young stars can produce CO shells 
around H II regions inside of a larger cloud or complex. 
This is a clear case of a small-scale offset between Ha 
and CO. However, the physical scales in Figure [3] are too 
large for this effect to have much impact, it should occur 
on scales more like ~10 pc. 

Similarly, the scales over which Td cp diverges between 
CO and Ha peaks (75 — 300 pc) are probably too large to 
be produced by drift between young stars and their par- 
ent cloud. A typical inter nal GMC velocity disper sion in 
M33 is a few km s _1 (ler; Rosolowsky et al. 2003). Over 
an average cloud lifet ime (~30 Myr; [Blitz_eLal] [20071 
IKawamura et al] I2009T) . this implies a drift of at most 
100 pc. This extreme case is just large enough to regis- 
ter in our plot but unlikel y to drive the signal we see at 
scales of 150 — 300 pc. See lEngargiola et all (|2003l) for a 
similar consideration of GMCs and H I filaments. 



Instead of drifts or offsets, what we observe is sim- 
ply a lack of direct correspondence between the CO and 
Ha luminosities of individual star-forming regions. The 
brightest 150 CO peaks are simply not identical to the 
brightest 150 Ha peaks. The bottom right panel of Fig- 
ure Q] shows this clearly; about a third of the peaks in 
M33 are nearer to another peak of their own type (i.e., 
CO to CO or Ha to Ha) than to a peak of the other type. 
Thus Figure [3] shows that the ratio of CO to Ha emission 
varies dramatically among star-forming regions. In this 
case the size scale on the £-axis in Figure [3] is actually 
a proxy for the number of regions inside the aperture. 
In M33, apertures of 75 pc diameter usually contain a 
single peak. At 150 pc, this is still the case ~70% of the 
time, and at 300 pc only a few regions are included in 
each aperture. 

Why does the ratio of CO-to-Ha vary so strongly 
from region-to-region? The efficiency with which gas 
form stars may vary systematically from region to region 
(with high Ha peaks being high-efficiency regions), star- 
forming regions may undergo dramatic changes in their 
properties as they evolve (with Ha peaks being evolved 
regions) , or the mapping of observables to physical quan- 
tities (Equations[T]and[2|) may vary from region to region. 

It is difficult to rule out region-to-region efficiency 
variati ons, but there is also no strong evidence for 
them. iLerov et al.l (|2008f) looked for systematic varia- 
tions in Tdcp as a function of a number of environmen- 
tal facto rs and found little ev i dence f or any systematic 
trends . IKrumholz fc McKed (|2005l) : IKrumholz fc Tan! 
(|2007|) suggested that the clou d free-fall time determine s 
Tdep to first order, but based on lRosolowskv et all (2003), 
the dynamic range in free-fall times for M 33 clouds is low. 
On the other hand, iGardan et all (|2007l ) found unusually 
low values of Td cp in the outer disk of M33. 

There is str ong evidence for evolution of star - 
forming regions iFukui et all (I1999D. iBlitz et all (|2007l) . 
IKawamura et all (|2009j), and iChen et al.l (|2010l ) showed 
that in the LMC, the amount of Ha and young stars 
associated with a GMC evolves significantly across its 
lifetime. In our opinion this is the most likely explana- 
tion for the behavior in Figure [3J Star-forming regions 
undergo a very strong evolution from quiescent cloud, to 
cloud being destroyed by H II region, to exposed cluster 
or association. When an aperture contains only a few 
regions, Td op for that aperture will be set by the evolu- 
tionary state of the regions inside it. That state will in 
turn determine whether the aperture is identified as a CO 
peak or an Ha peak. CO peaks will preferentially select 
sites of heavily embedded or future star formation while 
Ha peaks are relatively old regions that formed massive 
stars a few Myr ago. 

Region-to-region variations in the mapping of observ- 
ables (CO and Ha) to physical quantities (H 2 mass and 
SFR) are expected. Let us assume for the moment that 
the ratio of H2 to SFR is constant and independent of 
scale. Then to explain the strong scale dependence of 
the ratio of CO to Ha in Figure [3] there would need 
to be much more H2 per unit CO near Ha peaks and 
many more recently formed stars per ionizing photon 
near the CO peaks. At least som e of these effects have 
been claimed: e.g., Ilsraell (|1997l ) fi nd a strong depen- 
dence of Xco on radiation field and iVerlev et all (2010) 
suggest that incomplete sampling of the IMF in regions 
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with low SFRs drive the differences they observe between 
star formation tracers. However, both claims are contro- 
versial and it seems very contrived to invoke a scenario 
where only this effect drives the breakdown in Figure [3l 
It seems more plausible that the mapping of observables 
to physical quantities represents a secondary source of 
scatter correlated with the evolutionary state of a region 
(e.g., the age of the stellar population). 

5.1. Comparison to a Simple Model 

We argue that the behavior seen in Figures |3] and @] 
comes from averaging together regions in different states. 
Here we implement a simple model to demonstrate that 
such an effect can reproduce the observed behavior. 

The model is as follows: we consider a population of 
regions. We randomly assign each region to be an "Ha 
peak" or a "CO peak" with equal chance of each. CO 
peaks have 5 times as much CO as Ha and Ha peaks have 
5 times as much Ha as CO (roughly driven by the differ- 
ence between the results for 75 pc apertures in Table [1]). 
Physically, the idea is simply to build a population of re- 
gions that is an equal mix "young" (high CO-to-Ha) and 
"old" (low CO-to-Ha) . Dropping an aperture to contain 
only a young (CO peak) or old (Ha peak) region will 
recover our results at 75 pc scales by construction. Next 
we average each of our original region with another, new 
region (again randomly determined to be either a CO or 
Ha peak). We add the CO and Ha emission of the two 
region together, record the results. We then add a third 
region (again randomly young or old), and so on. 

The result is a prediction for the ratio of CO to Ha 
as a function of two quantities: 1) the number of regions 
added together and 2) the type of the first region (CO 
or Ha peak). Using the average number of regions per 
aperture listed in Table [1] and normalizing to an average 
depletion time of 1 Gyr, we then have a prediction for 
T^p as a function of scale. This appears as the diamond 
symbols and dashed lines in Figure [3l 

Given the simplicity of the model, the agreement be- 
tween observations and model in Figure [3] is good. Our 
observations can apparently be explained largely as the 
result of averaging together star-forming regions in dis- 
tinct evolutionary states. At scales where a single region 
dominates, the observed Td ep is a function of the state 
of that region. As more regions are included, Td ep just 
approaches the median value for the system. 

5.2. Tdep vs. Scale at Different Radii 

The star formation law apparently breaks down (or at 
least includes a large amount of scatter) on scales where 
one resolution element corresponds to an individual star- 
forming region. The spatial resolution at which this oc- 
curs will vary from system to system according to the 
space density of star-forming regions in the system. 

The surface dens ities of star forma tion and H2 vary 
with radius in M33 ()Hever et al.ll2004l) . This allows us to 
break the galaxy into two regions, a high surface density 
inner part (r ga i = — 2.2 kpc) and a low surface density 
outer part ( r ga i = 2.2 — 4.5 kpc). We measure the scale 
dependence of Td ep for each region in the same way that 
we did for all data. An important caveat is that the DIG 
subtraction becomes problematic for the outer region, 
removing a number of apparently real but low-brightness 



H II regions from the map. We achieve the best results 
for large radii with the DIG subtraction is turned off and 
report those numbers here. The basic result of a larger- 
scale of divergence in the outer disk remains the same 
with the DIG subtraction on or off. 

We find the expected result, that Td ep for CO and Ha 
peaks diverges at larger spatial scales in the outer disk 
than the inner disk. In both cases the ratio of Td ep at CO 
peaks to Td cp at Ha peaks is ~1 for 1200 pc apertures. 
For 600 pc apertures that ratio remains ~1 in the inner 
disk but climbs to —2 in the outer disk, suggesting that 
by this time there is already some breakdown in the SFR- 
H2 relation. For 300 pc apertures, the same ratio is —1.7 
in the inner disk and ~3 in the outer disk. It thus appears 
that at large radii in M33 the star formation law breaks 
down on scales about twice that of the inner disk, though 
the need to treat the DIG inhomogeneously means that 
this comparison should not be overinterpreted. 

6. CONCLUSIONS 

Our main conclusion is that the molecula r star forma- 
tion law observed in M33 at large scales (e.g. JHever et al.1 
120041 : IVerlev et af1l2010[ ) shows substantial scale depen- 
dence if one focuses on either CO or Ha peaks. The 
median depletion time (or CO-to-Ha ratio) measured in 
a 75 pc diameter aperture (derived from averaging —150 
such apertures) varies by more than an order of magni- 
tude between CO and Ha peaks. At large (—kpc) scales 
this difference mostly vanishes. We argue that the scale 
for the breakdown is set by the spatial separation of high- 
mass star-forming regions, with the breakdown occurring 
when an aperture includes only a few such regions in 
specific evolutionary states (a scale that corresponds to 
-300 pc in M33). 

In this case the scaling relation between gas and star 
formation rate surface density observed at large scales 
does not have its direct origin in an instantaneous cloud- 
scale relation. Individual GMCs and H II regions will 
exhibit a CO-to-Ha ratio that depends on their evolu- 
tionary state (likely with significant additional stochas- 
ticity) and as a result the —150 brightest objects at a 
given wavelength will be a function of what evolutionary 
state that observation probes. This div ergence is consis- 
tent w ith recent results from the LMC ([Kawamura et all 
l2009i ) indicating that individual GMCs exhibit a range 
of evolutionary states over their 20 — 30 Myr lifetime. 

This does not mean that comparisons of tracers of re- 
cent and future star formation on small scales are use- 
less. To the contrary, such observations contain critical 
information about the evolution of individual clouds as a 
function of time and location that is washed out at large 
scales (> 300 pc in M33). However, once one moves into 
the regime where a single object contributes heavily to 
each measurement, it is critical to interpret the results 
in light of the evolution of individual clouds. 
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APPENDIX 

Here we test how our method of selecting peaks and the removal of a diffuse ionized component affect our results. 

First, we repeating the analysis on the original maps without any DIG subtraction. We show the results in the left 
panel of Figure El along with the original measurements (in gray) from Figure [H To hrst order, the scale-dependence 
of Tdep is unchanged, but the H2 depletion are offset; Td ep derived from maps with no DIG subtraction is a factor of 
~2 smaller from the DIG-subtracted maps. 

As a second test we assess the impact of the particular choice of aperture positions. In main text we used a 
peak-finding algorithm. Here we test the effect of using published positions of GMCs and H II regions instead. The 
right panel in Figure [5| sho ws our origin a l dat a in gray while the red and blue points are Td ep derived using the 
iRosolowsky et al.r(|2007l ) and lHodge et al.l (|2002f ) catalogs. The median Td ep at large scale is unchanged from Figure [3] 
(gray values). However, for apertures centered on GMCs, Td ep on small scales does change from our analysis. This 
difference originates in different numbers and locations of the positions that are studied. In our original analysis, we 
study 172 positions which have CO emission peaks above 3c The Rosolowsky catalog, on the other hand, consists of 
only 140 positions inside a galactocentric radius of 4.5 kpc. In addition, a subset of the two samples targets different 
regions in M33: First, the catalog positions tend to be more clustered than the "peak" positions which leads to a 
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Fig. 5. — Tests of systcmatics. We repeat the measurement of T^ cp and remake Figure \3\ varying our methodology. Left: Results if we 
perfor m no DIG-subtraction o n the Ha maps. Right: Results if we use the known positions of GMCs and H II regions from Rosolowsk v et all 
(2007) and Hodge ct al. (2002). In both panels, our original results using the DIG-subtracted maps and an automated peak finder appear as 
gray points. Our qualitative conclusions are largely robust to these changes in methodology, though the overall normalization of CO-to-Ha 
ratio does depend on the DIG subtraction. 



somewhat larger number (5 — 15%) of objects in the smaller apertures and a smaller deviation in depletion times for 
CO or H II centered apertures. Second, while the molecular gas surface densities at the positions of the two samples 
do not differ significantly, the star formation rate surface densities are a factor of ~3 higher for the catalog positions 
as compared to the (more numerous) "peak" positions. This leads to shorter H2 depletion times on small scales for 
the catalog sample. 

Both tests show that the analyzed scale dependence of the star formation relation and the determination of its 
origin is not strongly dependent on the particular methodology chosen in this paper. While global shifts in the derived 
depletion times can arise due to the subtraction of diffuse Ha emission, we find only small variations in the scale 
dependence due to different selection of positions where we perform our study. 



